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ABSTRACT 


Extrapolations  of  numerical  data  obtained  from  finite  cluster  calcula¬ 
tions  to  the  infinite  volume  limit  can  often  give  incorrect  results.  We  discuss 
four  separate  cases:  (a)  the  intensity  of  the  lowest  two-photon  absorption  in 
the  infinite  polyene,  (b)  bond  alternation  in  the  infinite  polyene,  (c)  Cooper 


type  pairing  in  the  simple  Hubbard  model,  and  (d)  pairing  within  the  extended 
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I.  Introduction 

Approximate  analytic  techniques  have  been  known  to  lead  to  qualita¬ 
tively  incorrect  predictions  for  solid  state  systems  with  strong  electron  corre¬ 
lation.  Exact  numerical  techniques  are  therefore  often  employed  to  determine 
the  behavior  of  such  systems.  This  approach  consists  of  two  steps,  the  first 
of  which  involves  an  exact  calculation  for  a  finite  cluster.  The  second  step  in¬ 
volves  an  extrapolation  of  the  finite  cluster  data  to  the  thermodynamic  limit. 
We  point  out  that  extrapolations  are  not  always  valid,  and  in  certain  cases 
can  lead  to  erroneous  theoretical  predictions. 

We  believe  that  there  is  a  strong  necessity  for  such  a  discussion  cur¬ 
rently.  The  technique  of  using  exact  numerical  data  for  finite  systems  to 
determine  the  behavior  in  the  solid  began  with  the  classic  calculations  for  the 
Heisenberg  chain  by  Bonner  and  Fisher1 ,  who,  however,  emphasized  the  ne¬ 
cessity  of  having  one  dimension  and  short  range  interaction.  In  recent  years, 
however,  the  numerical  approach  has  been  employed  more  and  more  liberally, 
and  calculations  are  routinely  being  done  for  higher  dimensions,  as  well  as  for 
Hamiltonians  with  long  range  interactions.  Our  hope  is  that  the  present  work 
will  lead  to  more  systematic  studies  of  finite  size  effects. 

We  discuss  four  specific  cases  here.  The  correct  solutions  to  these 
problems  are  of  interest  to  chemists  and  physicists.  These  four  different  cases 
demJbstrate  different  finite  size  effects,  and  in  each  case  we  show  that  literal 
interpretation  of  finite  cluster  calculations  predict  qualitatively  incorrect  result. 

II.  Exact  Finite  Cluster  Calculations 

A,  The  intensity  of  the  twophoton  absorption  fTPAl  to  the  2Aff  and  other 
low  lying  Aa  states  in  the  infinite  polvene  The  occurrence  of  the  lowest  two¬ 
photon  state,  the  2A#,  below3  the  lowest  one-photon  state,  the  1  Bm  (note 
that  optically  inactive  Bn  states  need  not  be  considered),  in  finite  polyenes 
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necessarily  requires  a  correlated  description  for  the  Tr-electrons2,3 .  The  entire 
discussion  of  TPA  has.  however,  so  far  been  limited  to  the  energy  of  the  2 Ag 
state.  Here  we  extend  this  discussion  to  the  question  of  the  intensity  of  the 
TPA  to  the  2Ag  (and  other  low  lying  At  states)  as  a  function  of  the  chain 
length  N. 

This  discussion  will  be  in  the  context  of  the  extended  Hubbard  Hamil¬ 
tonian, 

H  =  +c+ ,  „c,,„) 

«,<f 

]CnUn.,i  1  -  1)  (1) 

i  i 

where  all  operators  and  parameters  have  their  usual  meanings.  We  begin  with 
the  U  =  V  =  0  limit.  In  Fig.  1  we  show  the  occupancies  of  the  highest  two 
valence  band  levels  and  lowest  two  conduction  band  levels  for  the  1  At  (the 
ground  state),  the  lfl»,  the  2At  and  the  2 Bu  for  arbitrary  N.  TPA  is  given 
by  the  third  order  nonlinearity  w;u\u>,  — u>)  =  Xrpx' 

TPA  iti  (wiB-  ~  w)(wMf  ~  2 u>)(u>{S.  -  w) 

where  fi  is  the  dipole  operator.  Conventional  wisdom  has  been  that  TPA 
occurs  every  time  an  energy  denominator  goes  to  zero  in  Eq.(2).  We  show 
below  that  this  is  not  true.  0 

The  intensity  of  the  TPA  to  the  2 At  was  calculated  exactly  for  U  = 
V  =  0  and  for  N  =  6,  8,  10  and  12,  typical  chain  lengths  for  which  calculations 
are  done  for  nonzero  U  and  V.  As  seen  in  Fig.  2.,  TPA  intensity  increases  with 
N  here,  suggesting  a  large  finite  TPA  intensity  for  N  — *  oo.  This,  however,  is 
incorrect. 

The  model  in  Eq.  (1)  has  charge-conjugation  symmetry,  which  implies 
that  for  every  single-particle  valence  band  level  at  energy  -e  there  exists  a 
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conductioD  band  level  at  e.  In  the  infinite  chain,  dipole  allowed  absorptions 
occur  only  between  these  symmetrically  placed  molecular  orbitals,  in  agree¬ 
ment  with  /^-conservation.  Thus  for  TPA  to  any  intraband  Ag  state,  only  two 
“symmetric”  £„  states  need  be  considered.  For  the  2Ag,  these  two  Bu  states 
are  the  l£u  and  the  2 Bu.  For  any  N,  the  absolute  values  of  the  transition 
dipole  moments  (2Ag\n\\B%)  and  (2Ag\fi\2Ba)  are  equal,  while  because  of  the 
fermion  character  of  the  electron  the  products  {lAg\n\lBa)(\Bu\ti\2Ag)  and 
(lAg\n\2Bu)(2Bu\n\2Ag)  are  of  opposite  signs.  At  N  — ♦  oo,  these  two  prod¬ 
ucts  become  equal  in  magnitude,  while  having  opposite  signs.  This  implies 
that  the  contribution  to  x^  Eq.(2)  by  the  2 Ag  vanishes4,  and  the  TPA 
intensity  to  the  2A^  and  other  band  edge  two-photon  states  is  zero  for  jV  — ►  oo. 

We  emphasize  that  the  above  analysis  does  not  imply  a  zero  intensity 
for  the  total  TPA  to  intraband  two-photon  states.  Even  with  the  restriction 
to  “symmetric”  Ba  states,  there  can  be  very  weak  TPA  due  to  relatively  high 
energy  Ag  states  that  are  slightly  removed  from  the  band  edge.  For  example, 
one  can  have  such  a  high  energy  intraband  At  state  that  is  coupled  to  the 
l£„  and  a  B%  state  that  requires  excitation  from  deep  inside  the  valence  band 
to  deep  inside  the  conduction  band.  Because  of  the  small  nonzero  energy 
difference  between  the  two  Ba  states  now,  the  cancellation  i3  no  longer  total. 
Weak  TPA  to  states  away  from  the  band  edge  is  expected,  and  explains  the 
shift  of  the  TPA  to  higher  energy  from  the  optical  band  edge  in  the  long 
chain3.  Nevertheless,  at  infinite  N,  TPA  to  band  edge  two-photon  states  goes 
to  zero,  in  contradiction  to  the  prediction  form  Fig.  2.  Since  calculations  of 
correlated  chains  are  limited  to  the  region  where  TPA  increases  with  N,  the 
integration  of  these  two  contradictory  results  is  an  important  issue. 

For  the  noninteracting  case  TPA  can  be  calculated  for  arbitrary  N. 
In  Fig.  3  we  show  the  TPA  intensity  to  2A#  as  a  function  of  N.  It  is  seen 
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that  the  TPA  intensity  does  go  to  zero  at  very  large  N,  but  only  after  a 
maximum  is  reached.  Thus  calculations  near  N  =  10  predict  qualitatively 
incorrect  results.  On  the  other  hand,  once  the  mechanism  is  understood,  one 
can  probe  the  N-dependence  of  |S~/S+|,  where  S~  and  S+  are  the  negative 
and  positive  contributions  to  the  TPA.  This  is  shown  in  Fig.  4,  where  the 
expected  N-dependence  is  observed. 

We  are  now  in  a  position  to  calculate  TPA  to  2 A,  in  long  correlated 
chains.  We  have  adopted  the  novel  configuration  interaction  (Cl)  approach 
of  Srinivasan  and  Ramasesha6  for  the  interacting  Hamiltonian  of  Eq.  (1), 
with  a  cutoff  4t  in  the  energies  of  single- particle  configurations.  The  approach 
can  give  the  2A9  below  the  1  Bu  even  in  long  chains.  In  Table  1  we  compare 
calculated  |5“/5+|  for  many  different  U,  V  and  6.  The  ratio  cannot  be  1  in 
N  =  20  -  30,  but  comparisons  with  the  noninteracting  cases  strongly  suggest 
that  the  TPA  to  the  2 Ag  should  be  very  weak  in  the  long  chain  limit  even 
for  nonzero  Coulomb  interaction,  in  contrast  to  what  would  be  predicted  from 
exact  evaluation7  of  the  TPA  near  N  =  10.  This  prediction  is  in  agreement 
with  the  tiny  TPA  to  2 A,  in  poly-BCMU  polydiacetylene8,  and  the  noted 
absence  of  the  two-photon  resonance  in  ^-carotene9. 

B.  The  bond-alternation  problem  in  the  one  dimensional  half-filled  band 
Within  the  Su-Schrieffer-Heeger  Hamiltonian, 

H  =  “  a(yi  ~  yi+>)](ckc'+».®  +  c*+i,<rc«>)  +  \KY^(y'  ~  y,+  1^  ^ 

i,e  i 

bond-alternation  in  the  infinite  polyene  is  unconditional.  Coulomb  interac¬ 
tions  can  enhance  this  bond  alternation.  For  long  range  Coloumb  interactions 
H<t  of  the  form  , 

Htt=uY<  "U".,l  +  Y,  *>("•  -  !)(".+;  -  1)  <4) 
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two  sharp  inequalities  give  the  ground  state  broken  symmetry10.  For  otf- 
Coulomb  parameters  that  obey  the  convexity  condition,  >  2VJ,th<- 

ground  state  can  only  be  a  bond-order  wave  (BOW)  with  alternating  bonds 
or  a  charge  density  wave  (CDW)  with  a  periodic  modulation  of  the  on-site 
charge  density.  The  dominant  broken  symmetry  is  given  by  the  inequality, 

X>„+,  £  +  (5) 

where  a  smaller  left  hand  side  (LHS)  gives  BOW,  while  a  smaller  right  hand 
side  (RHS)  gives'a  CDW.  All  P&riser-Parr-Pople  (PPP)  parametrizations  for 
the  infinite  polyene  then  predict  enhanced  BOW10. 

The  approach  involving  extrapolation  technique  was  also  employed  in 
this  case,  by  Ramasesha  and  Soos11.  These  authors  calculated  the  electronic 
energy  gained  for  fixed  bond  alternation,  when  compared  to  the  equal  bond 
length  case,  for  open  even  chains,  closed  N  =  4n  rings,  and  closed  N  =  4n  4-  2 
rings,  and  extrapolated  the  energy  gains  in  the  three  series  against  l/N  for 
the  PPP-Ohno  Hamiltonian.  Enhanced  energy  gains  compared  to  the  simple 
Huckel  case  was  taken  to  be  the  proof  of  enhanced  BOW.  While  the  result 
for  the  carbon-based  polymers  agrees  with  the  prediction  from  Eq.  (5),  let  us 
examine  this  approach  in  detail. 

The  extrapolation  of  the  4n-ring  series  is  invalid,  since  the  RHS  of  Eq. 
(5)  in  this  case  always  has  an  extra  term  that  precludes  the  CDW  entirely  and 
guarantees  enhanced  BOW.  The  even  chain  extrapolation  has  many  different 
finite  size  effects.  First,  even  chains  have  an  odd  number  of  bonds,  thus 
precluding  equal  bond  lengths  for  small  N.  Second,  the  range  of  the  Coulomb 
interaction  within  PPP  models  is  longer  than  the  N  for  which  calculations 
are  being  done.  This  leads  to  a  different  Hamiltonian  for  each  N.  Finally, 
real  space  many-electroo  configurations  have  entirely  different  matrix  elements 


6 


of  Hee  within  the  open  chain  boundary  condition  and  for  the  periodic  ring. 
For  example,  the  configuration  ...20202020...,  where  the  numbers  denote  site- 
occupancies,  has  a  matrix  element  4 U  —  8Vi  +  8V2  -  8F3  +  4Vt  in  the  ring, 
but  the  matrix  element  is  41/  —  4-  6V2  —  5V3  -f  4V4  —  3V5  4-  2V6  —  V7  in 

the  open  chain.  This  last  difference  shifts  the  boundary  between  the  BOW 
and  the  CDW  considerably,  such  that  the  open  chain  would  predict  enhanced 
BOW  even  as  Eq.  (5)  predicts  a  CDW.  Since  the  4n  +  2  ring  series  can  only 
have  a  few  data  points,  this  makes  the  extrapolation  technique  itself  suspect. 

We  demonstrate  the  above  by  doing  the  open  chain  calculation  for  the 
PPP-Ohno  parameters  corresponding  to  polysilane.  It  has  been  assumed  that 
this  has  a  simple  BOW  ground  state  with  very  large  alternation12.  We  choose 
a  U  of  9.04  eV  and  evaluate  Vj  from  the  PPP-Ohno  relationship, 


r  V 

3  "  (1  +  0.39J?2]1/2  (6) 

where  the  R}  are  the  distances  between  the  polysilane  “sites”,  and  are  taken 
from  literature12.  Eq.  (5)  predicts  a  CDW,  with  different  charge  densities  on 
the  orbitals  on  the  same  silicon  atom.  The  extrapolation  technique  neverthe¬ 
less  predicts  enhanced  BOW,  as  seen  in  Fig  5. 

Our  motivation  here  was  to  demonstrate  the  inapplicability  of  the  ex¬ 
trapolation  technique  for  the  bond  alternation  problem  and  not  to  claim  that 
polysilane  is  a  CDW.  Whether  or  not  polysilane  is  a  CDW  depends  on  whether 
the  PPP-Ohno  model  is  truly  appropriate  for  polysilanes.  It  is,  however,  inter¬ 
esting  to  note  that  the  occurrence  of  the  lowest  two-photon  absorption  consid¬ 
erably  above  the  lowest  one-photon  absorption13  occurs  naturally  within  the 
CDW  model. 

C.  Pairing  in  Hubbard  models  A  great  deal  of  numerical  effort  has  gone  into 
investigating  whether  pairing  of  the  Cooper  type  occurs  in  the  weakly  doped 
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two-dimensional  Mott-Hubbard  insulator.  More  recently,  such  a  pairing  mech¬ 
anism  has  been  proposed14  for  doped  C60.  The  two  dimensional  calculations 
have  now  been  extended  to  the  case  of  weakly  coupled  layers15.  Specifically, 
these  calculations  involve  determination  of  the  ground  state  energies  of  the 
half-filled  band  system,  and  the  electronic  energies  needed  to  add  one  and  two 
electrons  (holes).  Pairing  is  supposed  to  occur  if  the  quantity 

AE  =  2E{N  +  1)  -  E(N)  -  E{N  +  2)  (7) 

is  positive.  In  spite  of  considerable  effort,  it  is  still  not  entirely  clear  whether 
pairing  occurs  in  two  dimension  or  not,  or  what  the  actual  mechanism  of 
pairing  is  in  the  various  cases. 

We  focus  here  on  the  mechanism  of  the  pairing,  and  what  this  says 
about  pairing  in  the  infinite  system.  Fye  et  al1®  showed  that  pairing  of  the 
above  kind  is  not  limited  to  two  or  three  dimension,  but  can  also  occur  in  one 
dimensional  periodic  rings  with  N  —  An  atoms.  The  limitation  to  4n  rings  is 
very  suggestive,  as  small  4n-rings  undergo  Jahn-Teller  distortions.  In  Fig.  6 
we  show  our  results  of  the  pairing  calculations  for  N  =  4  and  8.  The  pairing 
disappears  for  the  Jahn-Teller  distorted  rings.  To  understand  this  better,  we 
calculated  the  effect  of  the  Hubbard  U  on  the  Jahn-Teller  distortion  of  the  4n 
and  4n  —  1  electron  systems.  These  restilts  are  shown  in  Fig.  7.  It  is  seen  that 
U  strongly  reduces  the  gain  in  energy  for  the  Jahn-Teller  distortion  of  the  even 
(4n)  electron  system,  but  the  distortion  of  the  odd  numbered  (4n  -1)  electron 
system  remains  virtually  unaffected.  This  then  gives  us  the  mechanism  of 
the  pairing.  The  information  regarding  the  tendency  to  distort  is  built  into 
the  wave  function  in  the  form  of  the  various  bond  order  parameters.  By  not 
letting  the  system  distort,  we  are  overestimating  the  energies  in  each  case,  but 
as  Fig.  7  indicates,  this  overestimation  is  much  stronger  for  the  odd  electron 
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case  than  for  the  even  electron  case.  The  calculated  energy  required  to  add 
one  hole  (or  electron  1  to  the  undistorted  system  is  therefore  too  high,  thereby 
giving  “pairin*  ’ 

This  same  mechanism  can  explain  “pairing”  in  higher  dimension.  In 
Fig.  8  we  show  the  results  of  our  calculation  for  the  simple  cube,  and  for  var¬ 
ious  distorted  forms.  For  the  undistorted  simple  cube,  we  merely  reproduce 
the  results  of  White  et  al.17,  —  two  distinct  regions  are  found  in  which  pairing 
is  seen  to  occur.  For  the  cube  stretched  along  one  axis,  we  see  that  the  pairing 
still  exists.  This  is  not  the  Jahn-Teller  mode  for  the  cube,  the  noninteracting 
limit  still  having  degeneracies.  However,  the  pairing  at  small  U  has  disap¬ 
peared  for  both  the  Jahn-Teller  distorted  cube  as  well  as  the  structure  with 
all  three  hopping  integrals  different.  This  clearly  indicates  the  relationship  be¬ 
tween  pairing  and  the  degeneracies  in  the  single  particle  limit.  We  have  once 
again  calculated  the  effect  of  the  correlation  parameter  U  on  the  Jahn-Teller 
distortions  .  The  result  is  identical  to  that  in  Fig.  7,  the  tendency  to  distort 
is  almost  unaffected  for  the  7  electron  system,  but  for  the  6  electron  system 
this  is  strongly  suppressed.  The  mechanism  of  pairing  at  small  U  in  the  cube 
is  then  the  same  as  in  the  4n  ring.  In  the  absence  of  distortion,  the  energy  to 
add  one  hole  or  electron  is  being  overestimated. 

We  have  repeated  such  calculations  for  a  number  of  three  dimensional 
molecules,  and  have  found  the  following.  For  an  N  site  system,  pairing  at 
small  U  occurs  only  if  the  ( N  —  1)  (or  (N  +  1))  electron  system  has  a  strong 
tendency  to  have  a  Jahn-Teller  distortion.  Thus  the  mechanism  of  pairing  at 
small  U  in  all  these  case  is  related  to  the  suppression  of  Jahn-Teller  distortion. 
Since  the  Jahn-Teller  distortion  results  from  the  discrete  level  degeneracies  in 
finite  molecules,  we  conclude  that  the  observed  pairing  in  in  the  small  U  region 
in  similar  calculations  is  a  finite  size  effect. 
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Fig.  8  explains  the  oscillation  in  the  pairing  energy  of  the  cube.  Pan  inn 
at  large  U  is  unrelated  to  the  tendency  to  have  Jalm-Teller  distortion.  We  have 
occasionally  found  such  pairing  at  large  U  also  in  other  systems,  though  not 
necessarily  accompanied  by  pairing  at  small  U.  The  occurrence  of  such  large  U 
pairing  in  nonbipartite  systems  (for  e.g.,  in  a  pentagonal  prism)  indicates  tha* 
this  is  also  unrelated  to  antiferromagnetism.  Currently  we  are  investigating 
whether  the  large  U  pairing  is  also  a  finite  size  effect,  and  if  so,  the  nature  of 
the  finite  size  effect.  The  present  work  indicates  the  absence  of  pairing  in  the 
two  dimensional  Hubbard  model,  as  well  as  in  weakly  coupled  layers. 

D.  Pairing  due  to  intersite  interactions  Many  calculations  exist  that  find  a 
pairing  driven  by  intersite  interactions  rather  than  the  on-site  Coulomb  inter¬ 
action.  Since  the  intersite  interactions  in  these  calculations  are  different,  it  is 
impossible  to  give  a  general  critique.  Nevertheless,  we  believe  numerical  cal¬ 
culations  here  would  in  most  cases  find  phase  segregation,  rather  than  simple 
pairing,  if  adequate  numbers  of  electrons  or  holes  are  added  to  the  half-filled 
band.  In  the  case  of  the  t-J  model,  this  tendency  to  phase  segregate  has  been 
discussed18.  We  choose  one  simple  example,  a  calculation  claiming  pairing 
driven  by  the  nearest  neighbor  Coulomb  interaction19  in  the  two-dimensional 
quarter- filled  band,  to  illustrate  our  point. 

For  large  U  and  V,  the  quarter-filled  band  is  a  CDW.  The  configuration 
with  one  extra  particle  is  as  shown  in  Fig.  9.  This  has  configuration  interaction 
with  a  higher  energy  configuration  which  is  at  energy  2V  relative  to  it.  If, 
however,  one  adds  two  particles  which  occupy  neighboring  empty  sites  in  the 
CDW,  the  lowest  energy  configuration  has  configuration  interaction  with  a 
configuration  that  is  at  an  energy  V  relative  to  it.  Thus,  Cl  should  give  pairing 
in  two  dimension,  though  not  in  one  dimension.  This  was  actually  numerically 
demonstrated  by  Mazumdar  and  Ramasesha19.  The  fallacy  in  this  argument 
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becomes  clear  once  three  particles  are  added  to  the  system  as  shown  in  Fig.  9 
Now,  the  lowest  energy  configuration  has  Cl  with  a  configuration  at  the  same 
energy,  so  that  a  cluster  of  three  particles  is  favored  over  a  cluster  of  two  and 
an  isolated  particle.  Similar  qualitative  arguments  favor  phase  segregation 
in  related  models.  Thus,  once  again,  thus,  the  pairing  found  within  such 
calculations  is  a  finite  size  effect. 

ID.  Conclusions 

The  major  conclusion  of  this  presentation  is  that  there  are  many  dif¬ 
ferent  finite  size  effects,  and  even  when  a  given  approach  is  satisfactory  for 
the  determination  of  one  property  of  a  given  class  of  materials  (e.g.,  the  en¬ 
ergy  of  the  2 Aa  in  the  infinite  conjugated  polyene),  it  does  not  necessarily 
mean  that  the  same  approach  is  satisfactory  for  other  properties  (in  this  case 
ground  state  bond  alternation  and  the  intensities  of  TPA).  Finite  size  effects 
should  be  investigated  for  each  different  problem  separately.  Probably  the 
most  useful  approach  to  avoid  the  pitfalls  is  to  have  a  physical  and  intuitive 
mechanism  of  the  enect  in  question  in  mind.  However,  this  can  hardly  be 
prescribed  as  a  general  approach.  One  simple  way  to  investigate  finite  size 
effects  is  to  carefully  investigate  the  size  dependence  of  the  property  in  ques¬ 
tion  within  the  noninteracting  model,  for  which  calculations  can  be  done  for 
arbitrary  sizes.  Even  though  this  procedure  is  simple  (and  an  obvious  step), 
finite  correlated  clusters  are  often  compared  to  the  infinite  noninteracting  sys¬ 
tem,  simply  because  the  latter  is  understood.  Such  a  procedure  can,  and  often 
does,  erroneously  ascribe  effects  due  to  finite  size  to  correlations. 
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Table  Caption 

Table  1.  The  ratio  of  the  negative  and  positive  contributions  to  the  2Ag  TPA 
in  long  polyene  chains  within  the  extended  Hubbard  model. 

Figure  Caption 

Figure  1.  The  single-particle  orbital  electron  occupancies  of  the  1  Ag,  15U, 
2 Ag  and  the  2 Bu.  Only  the  highest  (lowest)  two  valence  (conduction)  band 
levels  are  shown.  All  other  valence  (conduction)  band  levels  are  filled  (empty). 

Figure  2.  The  intensity  of  the  2 Af  TPA  in  short  finite  chains  within  the  Hiickel 
model.  Here  Et  is  the  15,  —  1  Ag  gap.  Notice  that  the  intensity  increases  with 
N. 

Figure  3.  The  intensity  of  the  2 Ag  TPA  as  a  function  of  N  within  the  Huckei 
model.  TPA  goes  to  zero  at  N  — »  oo. 

Figure  4.  The  ratio  of  the  negative  and  positive  contributions  to  the  2 Ag  TPA 
as  a  function  of  N  within  the  Huckei  model. 

Figure  5.  The  electronic  energy  difference  per  site  between  dimerized  and  the 
equal  bond  length  even  chains  for  the  parameters  of  Eq.  (6).  The  extrapolated 
N  -*  oo  values  lies  above  the  value  for  the  dimerized  infinite  Huckei  ring. 

Figure  6.  Pair  binding  energies  (see  Eq.  (7))  for  the  (a)  N  =  4  and  (b)  N  =  8 
rings.  In  each  case  the  solid  curve  corresponds  to  the  undistorted  ring,  the 
dashed  curve  to  the  Jahn-Teller  distorted  ring. 

Figure  7.  The  gain  in  electronic  energy  on  J ahn-Teller  distortion  for  the  N  —  4 
ring,  as  a  function  of  Hubbard  U .  The  solid  and  dashed  curves  correspond  to 
the  4  and  3  (5)  electron  cases  respectively. 

Figure  8.  Pair  binding  energies  in  the  simple  cube  and  three  distorted  forms, 
(i) . ,  simple  cube;  (ii)  - ,  stretched  cube,  tx  =  t2  >  t3\  (iii)  — 
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_ t  Jahn- Teller  distorted  cube;  (iv) - ail  three  hopping  integrals 

different.  Notice  absence  of  pairing  at  small  U  for  cases  (iii)  and  (iv). 

Figure  9.  Quarter-filled  extended  Hubbard  CDW  with  (a)  one  extra  particle, 
(b)  two  extra  particles,  (c)  three  extra  particles  (crosses  are  particles,  dots 
empty  sites).  CT  excitation  in  (a)  costs  energy  2V,  in  (b)  V,  and  zero  in  (c). 
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